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Quantum phase transitions in spin-boson 
systems: dissipation and light phenomena 

Karyn Le Hur 

Departments of Physics and Applied Physics, Yale University, New Haven, 
Connecticut, 06520, USA 

Spin-boson models are essentially useful in the understanding of quantum 
optics, nuclear physics, quantum dissipation, and quantum computation. We 
discuss quantum phase transitions in various spin-boson Hamiltonians, com- 
pare, and contrast them. We summarize the theoretical concepts and results, 
open questions and implementations of those ideas in cold atomic and QED 
cavity systems will also be addressed. The chapter is organized as follows. 

First, a large collection of harmonic oscillators (bosons) can simulate dissi- 
pation in quantum mechanics. Through a two-level system coupled to a bath 
of bosons, we investigate in detail the concept of "dissipation-driven" quantum 
phase transition. Another section will be devoted to the effect of dissipation 
(the bath of bosons) on the critical exponents associated with a well-known 
phase transition such as the disordered-ordered transition in the Ising model. 
Second, a spin-boson model can also describe the light-atom interaction. In 
particular, the Dicke model describing an ensemble of two-state atoms inter- 
acting with a single quantized mode of the electromagnetic field is well-known 
to exhibit a zero-temperature phase transition at a critical value of the dipolc 
coupling strength. Finally, we theoretically study the superfluid-Mott tran- 
sition of polaritons in the Jaynes-Cummings lattice system which consists of 
an array of coupled optical cavities each containing a two-level atom. 



1.1 Dissipative Transitions for the two-state system 

A dissipative two-state system generally refers to a two-level system coupled 
to a bath of harmonic oscillators (large collection of bosons) [1]: 




(1.1) 
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with, 

here, cr^ and tr z are Pauli matrices, A is the tunneling amplitude between 
the states with a z = ±1, and h the bias (detuning).* Moreover, Xi, p i: 
rrii, and are the coordinate, momentum, mass, and frequency of the ith 
harmonic oscillator. Here, Cj denotes the strength of the coupling with the 
ith oscillator. Information about the bath is encapsulated in the spectral 
function [1] J (to) = -| — 1 S(lo — u>i). This spin-boson model is a variant 
of the Caldeira-Leggett model [2] where the quantum system is a spin [1, 3]. 

1.1.1 Ohmic case 

In the case of ohmic dissipation, J(uj) = r]ue~ ii '^ c . It is then convenient to 
introduce the dimcnsionlcss dissipation (friction) coefficient a such that r/ = 
2ira. In fact, the emergence of a quantum phase transition can be understood 
from a perturbation theory in A/h where h > 0. More precisely, the order 
parameter obeys (a z ) + 1 <~ O(A) where A — (h/uj c ) 2a [4]. For a > 1, 
one observes that the spin tends to be trapped meaning that (cr z ) ~ —1 for 
h <C cj c . In contrast, for a < 1, A increases for smaller h and eventually 
reaches its maximum value .4. <~ 1. This indeed suggests the existence of a 
quantum phase transition at a c <~ 1. These impurity systems generally display 
both a classical (trapped) and quantum (untrapped) phase for the spin [5]. 

To better understand the nature of the dissipation-induced transition, one 
can integrate out the dissipative bath. This leads to an effective action which 
is reminiscent of the classical spin chains with long-range correlations [6] : 

S int = - f dr [ dT'a z {T)g(T-T')a z (T'), (1.3) 
Jo Jo 

with Q{t) oc 1/t 2 at long (imaginary) times uo^ 1 «t</? = l/ksT. It is rel- 
evant to observe that Anderson, Yuval, and Hamman [7] found an equivalent 
Ising model when studying the Kondo problem:^ 

H K = H kin + ^ ^2 (4:] c k'lS^ + 4i c fe'T» s ' + ) 

kk' 

+ Y E ( C I T ^'T - 4i c fe'l) S z + hS z , (1.4) 

kk' 



'Hereafter, the Planck constant will be set to h = 1 to avoid confusion with the bias h. 
t"We obtain the identity (S z ) = (o" z }/2. However, a similar relationship does not hold 
between (S x ) and (a x ). The electron operator must not be confused with the coupling 
constant Ci (or c^) in the spin-boson model which will prominently appear in the following. 
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where iJkin represents the kinetic energy of the electrons. This is maybe not 
so surprising: conduction electrons represent a dissipative bath and as a re- 
sult tunneling (spin flip) events are not independent. Instead, they feature 
long-range interactions in time and the equivalent Ising chain acquires long- 
range interactions. The trapped-untrapped transition in the spin-boson model 
is in fact equivalent to the ferromagnetic-antiferromagnetic transition in the 
anisotropic Kondo model [8, 9]. The equivalence between the anisotropic 
Kondo model and the spin-boson model with ohmic damping can be for- 
mulated rigorously through bosonization [10]. The untrapped region (for 
the spin) corresponds to the antiferromagnetic Kondo model J z > 0, while 
the trapped region corresponds to the ferromagnetic Kondo model J z < 
where the spin is fatally frozen in time. One gets the precise correspondence 
(pjjj — > A/lu c and (1 + 25/ir) 2 — > a, where p is the conduction electron 
density of states and the parameter 5 is related to the phase shift caused by 
the J z Kondo term and is given by S = tan _1 (— npJ z /A) [10]. Additionally, 
in the untrapped phase, the effective Kondo energy scale obeys [ll]:* 

Tk — A(A/D) a ^ 1 ~ a \ (1.5) 

for values of a not too close to the transition and close to the transition 
Tk assumes the exponential form of the isotropic, antiferromagnetic Kondo 
model; lnT^- oc l/(a c — a). The critical line separating the trapped and 
untrapped phase in the spin-boson model corresponds to a c = 1 + 0(A/lo c ). 
Next, we discuss recent developments on spin observables and dynamics. 



1.1.2 Exact Results 

As shown by Thouless, for the 1/t 2 Ising chain the magnetization (a z ) is not 
continuous at the transition [12]. This is also consistent with the fact that in 
the ohmic case, the phase transition is described by Renormalization Group 
equations similar to those in the XY model in two dimensions. Following 
Anderson and Yuval, the order parameter (a z ) jumps by a non-universal 
amount yfl/a c along the quantum critical line, a c = 1 + 0(A/lu c ) [13]. 

In the untrapped phase (a < a c ), one can apply the Bethe Ansatz approach 
to compute observables exactly [4, 14, 15]. For h <C Tk, we obtain [4, 15]: 

2 e 2<T^) ni + 1/(2 - 2a)] / h \ . „, 

= - ^ r[l + a/(2-2a)] fe) ' (L6) 

where b = alna + (1 — a) ln(l — a). Note that (a z ) oc H/Tk at small h, in 
keeping with the Kondo Fermi liquid ground state [1] . The local susceptibility 
of the spin converges to 1/A for a — > § in accordance with the two-level 



'We have introduced a high-energy cutoff D for the conduction electrons which is of the 
order of the Fermi energy; the relation between u) c and D is given, e.g., in Refs. [4, 14]. 
§For cv — ► 0, r[l] = 1, r[3/2] = i/tt/2, exp(6/(2(l - a)) = 1, and T K = A. 
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description and diverges in the vicinity of the phase transition due to the 
exponential suppression of Tk ■ Note that the longitudinal spin magnetization 
(the spin order parameter) (o~ z ) only depends on the "fixed point" properties, 
i.e, this is a universal function of h/Tx in the untrapped phase. Finally, in 
the trapped phase, one predicts (a z ) w — 1 + 0((A/oj c ) 2 ) [16]. 

The leading behavior of (a x ) in the untrapped phase is [4, 14, 15]: 

Urn (a x ) = -±-£- +M (a&, (1.7) 

h<s.T K la — 1 uj c A 

with 

e -6/(2-2a) r[l-l/(2-2q)] 

y/ir(l — a) 1 [1 — a/(2 — 2a)J 

As a — > 0, 7k — > A and M(0) = 1, so we check that (<r x ) — * 1. As we 
turn on the coupling to the environment, we introduce some uncertainty in 
the spin direction and (a x ) decreases. Note that (a x ) does not only depend on 
the fixed point properties; more precisely, in the untrapped (and highly non- 
perturbative) regime, (a x ) still contains a perturbative part in A/uj c stemming 
from the trapped region! For a < 1/2, the monotonic decrease of Tk/A 
dominates. In contrast, for a > 1/2, the first term in Eq. (1.7) dominates: 

(o'x)q>1/2,/ 1 ^0 = o 7 ■ (1-9) 

This result can also be recovered using a perturbation theory in A/ui c [4] . This 
emphasizes that the observable (cr x ) is continuous and small at the quantum 
phase transition. This is also consistent with the work by Anderson and Yuval 
which predicts (cr x ) ~ A/ui c exactly at the quantum phase transition [13]. 

Finally, we can also check that the spin component (o~ x ) evolves continuously 
close to a = 1/2. In the limit a — > 1/2, one can take M(a) — (4/-7r)r(l — 
2a) -> 4/(?r(l - 2a)) and use the identity D{a = 1/2) = 4lu c /it [4, 14] to find 



(a x ) - -(4/Tr)yjT K /Dln(T K /D), (1.10) 
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in agreement with the "non-interacting" resonant level description valid at 
the specific point a = 1/2 [4]. Exact results for the spin observables obtained 
using Bethe Ansatz are summarized in Fig. 1.1. Usually, a strong reduction of 
the off-diagonal diagonal elements of the spin reduced density matrix traduces 
quantum decoherence. Using the results above, we observe that quantum 
decoherence is prominent for a > 1/2 (where (a x ) becomes tiny and the 
entanglement between the spin and the bath becomes almost maximal [4, 15]). 



1.1.3 Spin dynamics and Entanglement 

Another useful quantity in the context of macroscopic quantum coherence is 
the occupation probability P(t) = {a z (t)}, where the system is subject to the 
non-equilibrium initial preparation o~ z (t = 0) = +1 and the initial density 
matrix is in a factorized form [1]. At time t = 0, the dynamics starts out. 

To study the spin dynamics it is convenient to perform a polaronic transfor- 
mation U = exp(— ia z fl/2) where f2 = ^^Ci/wiiUi^pi, such that the trans- 
formed Hamiltonian H' = U^HU takes the precise form [1] 

In the Heisenberg picture, the equations of motion for o~±(t) are easily ob- 
tained. Integrating and substituting them into the equation of motion for the 
transverse polarization a x (t), then one gets the exact formula: 

& z (t) = -^A 2 f (e- in ^e in ^a z (t') + a z (t') e - m ^ 'e in «) dt' . (1.12) 

On the other hand, to solve this equation, one usually uses approximations 
[17]. The first approximation generally consists to insert the free bath dy- 
namics when computing the commutator: 



[Cl(t)Mt')}=iYl [ sinMt-t'))- (1-13) 



3 



The next step is to average (1.12) with respect to the bath and to decouple 
the environmental exponentials from the spin. Using that: 



c 2 



(ft(i)fi(i') + Q(t')Q(t)) = J2 — 3 — coth ) cos(wj(i - t')), (1.14) 

j m j UJ j \2 / 

this leads to the evolution equation [17]: 

P(t)+ j F{t-t')P{t')dt' = 0, (1.15) 
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where the function T obeys T{t) = A 2 cos (Ai(t)) exp — (A 2 (t)), and 

AUt) = - / sm(wt)-^-du> (1.16) 

71" Jo ^ 
„ , ^ 1 f + °° , u , //M J(W) , 

^(i) = - / (1 - cos(wi)) coth — -^-du. 

K JO V 2 / w 

Through the Laplace transform one obtains (C denotes a Bromwich contour): 



c dX ' x, TTm- < L17 » 



At zero temperature and in the scaling limit A/uj c <C 1, one finds [1]: 



A R N 



T(X) = A e , (1.18) 

i 

where A e = A r (cos(7ra)r(l — 2a)) 2 < 1 -°) ; we have introduced the rcnormal- 
ized transverse field A r = A(A/w c ) a / 1 " Q! which is proportional to the Kondo 
energy scale Tk- This expression of P(t) coincides with the formula of P(t) 
obtained via the Non-Interacting Blip Approximation (NIB A) [1]. 

For a — > 0, one recovers perfect Rabi oscillations P(t) = cos(At) whereas 
for a — 1/2 one gets a pure relaxation P(t) = exp — (TrA 2 t/(2uj c )), which is 
in accordance with the non-interacting resonant level model [1]. For < a < 
1/2, the spin displays coherent oscillations (due to a pair of simple poles) 
leading to P co h(t) = acos((t + <p) exp(— jt) with the quality factor [1]: 

C cotf^V (1.19) 



7 \2(l-a) 

This quality factor has also been found using conformal field theory [18]. 

Recently, we have developed a time-dependent Numerical Renormalization 
Group (NRG) which allows us to confirm these results [19]; see Fig. 1.2. 
In particular, we have checked the expression given in (1.19) for the quality 
factor. In fact, the NRG represents a powerful theoretical tool to study those 
spin-boson systems [20, 21, 22]. Additionally, we have obtained the same 
quality factor using an exact analytical extension of the NIB A approach [23] . 

It is relevant to note that the coherent-incoherent crossover, corresponding 
either to the strong suppression of the off-diagonal elements of the spin re- 
duced density matrix (a x ) <~ A/u c — > or to the complete vanishing of the 
Rabi quantum oscillations, can be identified to the Toulouse limit a = 1/2 
and not to the quantum phase transition. It is also interesting to underline 
the correspondence between the prominence of spin-bath entanglement and 
the emergence of quantum decoherence [4, 15]. The entanglement entropy S 
between the spin and the environment displays a plateau at maximal entan- 
glement for 1/2 < a < 1. Using the time-dependent NRG [19], we have also 
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studied the strong coupling regime and in particular the crossover from inco- 
herent decay to localization at the quantum phase transition (Fig. 1.2). In 
particular, for 1/2 < a < 1, the spin dynamics remains purely incoherent. The 
authors of Ref. [24] report a multi-exponential form P(t) = J2n=i a n^ 2hnt ■ 
In contrast with the results of conformal field theory [18], the rate constants 
appearing in the multi-exponential fit, are not integer multiples of a single 
rate. Our NRG results rather support the form Pit) = exp[(— iTft-/2) a =] with 
the prerequisite that a\ii = 1. The exponent a a evolves linearly with a — 0.5. 

1.1.4 Sub-ohmic case 

Now, we focus on the sub-ohmic situation J(u>) = 2iraojl~ s oj s with < s < 1 
which exhibits a second-order quantum phase transition [5] by analogy to 
classical spin chains [6, 25]: in this case Q{t) oc 1/t 1+s in Eq. (1.3). The 
second-order quantum phase transition separates a localized (trapped) phase 
for the spin at small A from a delocalized (untrapped) phase at large A. 

For a second-order impurity quantum phase transition we can apply the 
following scaling ansatz for the impurity part of the free energy [5] , 

F mp = TFQA - Acl/T 1 /", hT- b ), (1.20) 

where we have re-introduced the detuning h and A c is the value of the trans- 
verse field at the quantum critical point. The critical exponent b should not 
be confused with the parameter b of Sec. 1.1.2. There is no independent dy- 
namical critical exponent for (0 + 1) dimensional models, formally z — 1. At 
zero temperature, the crossover from the quantum critical regime to one or 
other of the stable regimes, defines an energy scale h* which vanishes at A c , 
h* oc | A c — A\ hv . In a similar way, at finite temperature, we define the energy 
scale T* = |A C — A^. It should be noted that the ansatz (1.20) is usually well 
justified when the fixed point is interacting [5]; for a Gaussian fixed point the 
scaling function would also depend upon dangerously irrelevant variables. 
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Both analytical arguments, based on the equivalence to a 0(1) 4 theory 
and numerical simulations for the one-dimensional long-range Ising model 
show that the upper-critical dimension is d u = 2s [26, 27]. In other words, 
the transition obeys non-trivial critical behavior for 1/2 < s < 1 and the 
fixed point is interacting. More precisely, hyperscaling relations imply that 
there are only two independent exponents, e.g., v and b. A Ward identity 
for the spin-boson model ensures that b = (1 + s)/2 [4, 28, 29]. Defining the 
exponent 5 as (a z )(h,A c ) oc l/il 1 ^ 5 , one also finds 5 = (1 + s)/(l — s). This 
implies that the local susceptibility d(a z )/dh at the quantum critical point 
diverges as T~ s . Hyperscaling also guarantees an uj/T scaling, x" oc (A c ,T = 
0,w) cx |a>| _s sgn(w). Notably, this was found to be the exact decay exponent 
of the critical spin correlations in the long-range Ising model for all s [26]. 
Additionally, the exponent v diverges as l/y2(l — s) near the ohmic point s = 
1 [5]. For < s < 1/2, the transition is mean- field like and critical exponents 
obey, e.g., S = 3 and v = 1/s [26, 27]. Hyperscaling is violated. Very recently, 
using a powerful continuous time cluster Monte-Carlo algorithm, Winter et al. 
have shown that the quantum-to-classical mapping is valid in the case of the 
sub-ohmic spin-boson model [30]. On the other hand, note that the presence 
of a dangerously irrelevant variable for s < 1/2 impedes the correct extraction 
of the critical exponents with current versions of the NRG method. 

In Refs. [4, 28], we have studied in detail the entanglement properties in the 
sub-ohmic spin-boson model. In particular, for those second-order quantum 
phase transitions, the entanglement entropy exhibits a visible cusp at the 
quantum critical point and quantum decoherence becomes maximized. 

1.1.5 Realizations 

The spin-boson model can be realized in noisy charge qubits built of meso- 
scopic quantum dots or Cooper pair boxes [31, 32]. The gate voltage controls 
the detuning h and A corresponds to the tunneling amplitude between the dot 
and the lead(s) or the Josephson coupling energy of the junction. If the gate 
voltage source is placed in series with an external resistor, which can be mod- 
elled by a long LC transmission line, this may describe the spin-boson model 
with ohmic dissipation [14]. A one-dimensional Luttinger reservoir could also 
be used [33, 34]. The spin-boson model can also be derived when coupling a 
quantum dot to a boson and a fermion bath [35]. These nano-systems may 
also allow to address new important issues such as the non-equilibrium trans- 
port properties at a given quantum phase transition [36] . The sub-ohmic case 
s = 1/2 can be engineered through an RLC transmission line. Charge mea- 
surements could provide the quantity (cr z ), which represents the occupation 
of the dot or island. In a ring geometry, the application of a magnetic flux 
generates a persistent current which is proportional to (a x ) [14]. Solid-state 
two- level systems usually feature a coupling strength much below a = 1/2. A 
very promising candidate is the ultracold-atomic quantum dot coupled to a 
Bose-Einstein Condensate (BEC) [37], which allows an unprecedented control 
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of the coupling(s) between the qubit and the reservoir. The spin-boson model 
can also be engineered in trapped ions arranged in Coulomb crystals [38] . 



1.2 Dissipative spin array 

Now, we investigate the ground state of a spin array coupled to a common 
(large) collection of harmonic oscillators. We intend to show that when the 
coupling is longitudinal the system can be mapped onto a dissipative quantum 
Ising model. To simplify the discussion, we consider a one-dimensional chan- 
nel model for the bath and allow bosons (phonons, sound waves in a BEC, 
or photons) to propagate along a single direction with wavevector k and dis- 
persion ujk = vk. More general results are shown in Ref. [39]. Similar to the 
single spin case, the interaction between the spins and the boson bath reads, 

i k 

whereas the boson bath Hamiltonian reads Hb = ^2k u>ka k ak ( we se ^ ^=1). 
In this equation, Xi correspond to the positions of the spin impurities. 



1.2.1 Boson-mediated magnetic interaction 

First, we show that an exchange interaction between the spins is induced by 
the bosonic environment; this is analogous to the Rudcrman-Kittcl-Kasuya- 
Yosida interaction induced by a fermionic bath. For this purpose, we set 
the transverse field A = 0. In fact, the spin-boson interaction can be exactly 
eliminated through a polaronic (unitary) transformation along the lines of the 
single spin case [1]. More precisely, we perform a unitary transformation V = 

cxp^4 = cxp {^2n J^k 4 ke * fexi ( a fc ~ o}_ k )ai^j . The transformed Hamiltonian 

then takes the general form (where H = H Int + H B + J2i \ a iz)'- 

e- A He A = H+[H,A] + ^[[H,A},A} + ... . (1.22) 

Now, we choose the coefficients Ak = Ck/wk such that the induced term 
[i?B,^4] exactly cancels Hj nt . On the other hand, the transformed Hamilto- 
nian also produces an effective interaction between spins: 

5H = ~\ E E ^e ik ^~^a iz a jz . (1.23) 

1,3 k 
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If one envisions a cold atomic spin array coupled to a BEC reservoir, one can 
check that the factor c\/iVk is k- independent [39]. Furthermore, we identify: 

yyfcCx,-*,) = Jl^^ _ x ) /^ h ). (i.24) 

we have introduced the function sinc(x) = sin(a:)/x and £/, = w/w c (w c repre- 
sents the ultraviolet cutoff for the sound modes in the BEC and L denotes the 
length of the BEC). Notice that the induced interaction decays very rapidly 
for separations larger than the healing length ^ ■ In a realistic cold atom ex- 
periment, the distance between (atomic) quantum dots is of comparable size 
as the healing length. Therefore, we can restrict ourselves to nearest-neighbor 
spin interactions. It should also be noted that the effective interaction —K < 
is ferromagnetic and independent of the length L since c\juj k ocl/L [39]. 



1.2.2 Solvable dissipative model 

On the other hand, the boson bath is also expected to give rise to long-range 
correlations in time similar to the single spin case. In fact, to find those long- 
range correlations in time, it is judicious to integrate out the phonon (sound 
wave) modes using coherent state functional integrals. At a general level, this 
produces the following extra term in the action of the spin array [39]:' 

(1-25) 

where j3 — (fc^T) 1 , <Tfc(r) = °~jz(T)e lkX: > ', and 0„ = 2irn/[3 are Matsubara 
frequencies. Interestingly, the coupling to the boson bath provides two distinct 
contributions, which can be identified using the decomposition, ^(r)^^') = 
\{o- k {r)a* k {T) + a fe (r'K(r')) - \{a k {r) - a fc (r'))K(r) - a*(r')) [1]. 

The first term which is local in time is dominated by the Matsubara fre- 
quency term n — 0. This allows us to recover the ferromagnetic Ising contri- 
bution SH found above by applying the unitary transformation V. 

Now, similar to the case of a single two-level system coupled to a bosonic 
bath [1], the second (dissipative) contribution stems from the pole at uj k = 
—ifl n (or LOk = iSl n )- One can check that the main dissipative contribution 
is for i — j. This term produces an on-site long-range correlation in time 
oc 1/t 2 . Note that correlations in space are short-range, as a result of Eq. 
(1.24), whereas the correlations in time are long-range due to dissipation. 

Re- introducing the transverse field A, then we find that the spin dynamics 



use the convention cr fe instead of a kz to save space. 
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is dictated by the following action [39] : 



S, 



spin — 




(1.26) 



8 i-Jo 

dr a iz (T)o- jz (T). 



1 {' 

a5Z / d-rdr'air - T')a iz (T)o- lz (T') 



We have properly defined: 



a 



1 f°° 

(r — r') = — / dojJ{uj)e 



I" Jo 



(1.27) 



and the spectral function reads J(w) = ^SfeCfe^^ — w fe)- I n f act i one can 
establish that J(w) = 27ra<x> d where (i denotes the dimensionality of the BEC. 

Assuming we consider two-state systems close to the degeneracy point h = 
0, then one recognizes the action of a dissipative quantum Ising model [40]. 

1.2.3 Dissipative 4 theory 

For A = the system is in a ferromagnetic (ordered) phase whereas for 
large A » K the system will be in a paramagnetic (disordered) phase. The 
nondissipative quantum Ising model (a — 0) exhibits a second order phase 
transition around K ~ A, separating the paramagnetic and ferromagnetic 
phase, and can be mapped exactly onto the classical Ising model in d e ff = 
d+z = d+1 dimensions; z = 1 is the dynamic critical exponent that describes 
the relative dimensions of (imaginary) time and space. The classical model 
itself can be described by a <j) A theory. In particular, the quantum Ising chain 
lies in the universality class of the two-dimensional classical model [41]. 

When a > 0, the critical behaviour is most profoundly changed in the ohmic 
case (d = 1). First, as a reminiscence of the single two-state system, the bath 
will renormalizc the transverse field as A r = A (A/uj c ) a ^ 1 [1], such that 
the transition occurs at a smaller K <~ A r . Second, for d = 1, the dissipation 
will generate a term |u;|</>(g, iu) 2 in the effective </> 4 theory which will affect the 
critical exponents and change the universality class of the phase transition. 

1.2.4 Critical exponents 

Here, we summarize the main results for the critical exponents. In fact, the 
latter have been shown to be independent of the value of the parameter a 
[40, 42] and they have been thoroughly derived through a dissipative (f> 4 theory 
in d = (2 — e) dimension (e = 1) [42, 43] and Monte Carlo simulations [40]. 
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TABLE 1.1 



Critical exponents of (dissipative) quantum Ising models. 



Critical Exponents 


a 


7 


S 


V 


z 


Dissipative Quantum Chain 


0.319 


1.276 


~ 5 


0.638 


1.98 


d = 2 Quantum Ising model 


0.325 


1.241 


4.82 


0.630 


1 


Non-dissipativc Quantum Chain 


1/8 


7/4 


15 


1 


1 



For d = 1, due to the dissipative term \ijj\(j)(q,iuj) 2 in the effective 0(1) <p 4 
theory, the dynamic critical exponent is equal to z w 2, whereas z = 1 for 
d = (2,3). Therefore, the one-dimensional dissipative and two-dimensional 
quantum Ising models should behave similarly, since they have the same di- 
mension d e ff = d + z i=a 3. For d = 1, the </> 4 theory (up to second order in 
£•) predicts z « 1.98 [43] which is in accordance with the Monte Carlo simula- 
tions [40, 42]. A summary of some critical exponents for the one-dimensional 
dissipative (quantum) case and the two-dimensional quantum Ising model is 
presented in Table 1.1. The exponent 5 is defined in Sec. 1.1.4. Additionally, 
a iz cx | A — A c |^ and \ — da iz /dh cx |A — A c |~ 7 . The correlation length expo- 
nent obeys v sa 0.638 [40, 42]. In addition, from scaling laws, 7 = zv 1.276 
[40, 42, 43] and /3 = is/2& 0.319. One must also satisfy z = (6 - l)/2. 

1.2.5 Realizations 

In Ref. [39], we have shown that a spin-boson mixture of cold atoms can 
be used to engineer the quantum Ising model in a dissipative bath. We em- 
phasize that this setup embodies the first tunable realization of the quantum 
Ising model in a dissipative bath, and that several critical exponents can be 
measured using standard imaging techniques. On the other hand, the effect of 
a (nuclear spin) bath on the quantum phase transition of an Ising ferromagnct 
in a transverse field has been recently addressed experimentally in Ref. [44]. 

Regarding our spin-boson mixture [39] , the first specie lies in a deep optical 
lattice with tightly confining wells and forms a spin array; spin-up/down cor- 
responds to occupation by one/no atom at each site. The second specie forms 
a superfluid reservoir. Different species are coupled coherently via laser tran- 
sitions and collisions. Whereas the laser coupling mimics a transverse field for 
the spins, the coupling to the reservoir sound modes induces the ferromagnetic 
(Ising) coupling as well as dissipation. By measuring the critical exponents, 
one may confirm that the dissipative phase transition is still of second-order 
type where they are related by (3(5 — 1) = 7. The order parameter (cjj z ) now 
must go continuously to zero at a second-order phase transition which is in 
striking contrast to the case of a single two-level system coupled to a ohmic 
bath. Finally, for the dissipative spin chain, the value of the dynamic critical 
exponent can be directly obtained from the equality z = (6 — l)/2 ~ 2. 

Open issues in the field include the entanglement properties as well as the 
spin dynamics of the dissipative quantum Ising chain. Concerning the last 
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point, a progress in this direction has been recently achieved in Ref. [45]. 
However, the authors have neglected the long-range spin correlations in time 
induced by the bath assuming that the bosons have a nonzero inverse lifetime. 



1.3 One-mode superradiance model 

An ensemble of N two-level atoms interacting with a radiation field has been 
studied by many authors [46, 47, 48]. In the celebrated Dicke model, the atoms 
are assumed to be at fixed positions within a linear cavity of volume V and the 
separations between the atoms are large enough so that the interaction among 
them can be ignored. However, the fact that the atoms interact with the 
same radiation field, they cannot be treated as independent. The importance 
of treating the radiating atoms as a single quantum system was recognized 
by Dicke [46] . An exact solution for the Hamiltonian of N identical two-level 
atoms interacting with a single-mode quantized radiation field at resonance 
was given by Tavis and Cummings [47] . The thermodynamic properties of the 
system in the limit N, V — > +oo, N/V w finite have been first obtained by 
HcppandLieb [49]. They reveal a second-order classical phase transition. The 
single-mode Dicke model also admits a second-order quantum phase transition 
[50, 51]. Below, we analyze this second-order quantum phase transition. 

1.3.1 Hamiltonian 

We study the quantum regime of the one-mode superradiance (Dicke) model 
[46] where collective and coherent behavior of the pseudospins (atoms) is 
induced by coupling — with interaction A — to a physically distinct single- 
boson subsystem. In the following, when we refer to the Dicke model we shall 
mean the single-mode Hamiltonian unless otherwise stated. With omission of 
the A 2 term for the electromagnetic field, the Hamiltonian reads [51]: 

H = lu s J 7 + ua'a + -^= (a + a f ) (J+ + J J), (1.28) 

where this form follows from the introduction of a collective spin operator of 
length j = N/2. The resonance condition is w = w s . The thermodynamic 
limit of N — > +oo is thus equivalent to making the length of the pseudo- 
spin tend to infinity j — > oo. Here, uj s is the frequency splitting between 
the atomic levels, u is the frequency of the field mode, and k is the dipolc 
coupling strength. In addition, the collective atomic operators satisfy angular 
momentum commutation relations [J+, J_] = 2J Z and [J±, J z ] = I fJ±. 

The Dicke model is usually considered in the quantum optics approach of 
the rotating-wave approximation, which is valid for small values of A, and 
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consists to neglect the counter- rotating terms a' J+ and aJ_. This makes the 
model integrable and simplifies the analysis. Below, we follow Emary and 
Brandes, and derive exact results without the rotating-wave approximation 
[51]. In the thermodynamic limit (N,j) — > oo, the system shows a quantum 
phase transition at a critical coupling A c = y/ujuj s /2, and the system changes 
from a normal phase to a superradiant one. Superradiance means the decay 
of an excited population of atoms via spontaneous emission of photons. 

In fact, it is instructive to note that the problem reduces to a two- mode 
problem by using the Holstein-Primakoff transformation [52] of the angular 
momentum operators J z = (b^b — j), J + = b^ \j2j — Wb, and J_ = jt ; here, 
b and 6+ represent Bose operators [b,tf] = 1. The Hamiltonian becomes [51]: 



H = uj s (b^b-j)+iua^a + X(a + a 1i ) I b\ 1 - — + J I - —b) . (1.29) 




1.3.2 Normal phase 

In the normal phase (A < A c ), one can expand the square roots and this gives: 

H n = cj s b f b - juj s + wa)a + A(o + a r )(b + & f ). (1.30) 

This problem of two-coupled harmonic oscillators is exactly solvable. After 
diagonalizing the problem, one gets two independent (effective) oscillators: 

H n = E^cld + 4 n) <4c2 + \ (4 n) + E - n) -v-ws)- 3^s. (1.31) 



The bosonic operators (c\, c\, C2, c3,) are linear combinations of the original 

(n) 

operators and describe collective atom- field excitations. The energies E K ±' of 
the two independent oscillator modes take the form [51]: 

(e { ± i) ) 2 = \ (u 2 + lu 2 s ± y/(ul - W 2 ) 2 + 16A 2 CO^) . (1.32) 

It should be noted that the excitation energy E M is real only if lo 2 + us 2 s > 
yj (lo 2 — lo 2 ) 2 + 16\ 2 lulu s implying A < = A c . This underlines that 

H n is valid only for A < A c , i.e., in the normal phase. For large j, the ground 
state energy is —jui s whereas the excitation energies E± are 0(1). 

In the normal phase there is a conserved quantity (parity) which is given 
by: II = exp(27riV e ) and N e = a^a + J z + j represents the excitation number. 
It counts the total number of excitation quanta in the system and possesses 
two eigenvalues ±1. After the Holstein-Primakoff transformation, the parity 
operator becomes II = exp (i7r[a^a + b^b}) . In this formulation, there is an 
apparent analogy with the parity operator of a two-dimensional harmonic 
operator. The ground state has a positive parity with an even excitation 
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number: this is obvious at A = since the excitation number is zero. In 
addition, as the energy levels in this phase are non-degenerate the continuity 
of the ground state with increasing A ensures that it always has positive parity. 

1.3.3 Superradiant phase 

In the superradiant phase (A > A c ) the field and the atomic ensemble acquires 
macroscopic occupations, thus one has to redefine [51]: 

b ] ^S-s/n, (1.33) 

where m and n are 0{j). The problem can be still diagonalized, and again 
one gets a theory of two decoupled oscillators: 

H s = E^eW + E^e\e 2 - j + ^) (1-34) 

The expressions of the bosonic creation and annihilation operators (ei, e\, 62, e\) 
can be found in the Appendix A of Ref. [51] and the oscillator energies obey: 



where the parameter [i is given by ji = lolu s /4X 2 = X 2 /X 2 . Again, one can 

(s) 

check that the excitation E_ remains real as long as A > A c , and the ground 
state energy is given by — j[(2X 2 /uj) + (uj 2 ui /8X 2 )]. It is relevant to observe that 
each of every level of the total spectrum is doubly degenerate above the phase 
transition. In particular, one might also define a 1 * — > c 1 — y/m, — > $ + ^n, 
leading to the same energy spectrum. It implies that the symmetry of the 
ground state defined by the parity operator II is spontaneously brokem at 
A c . Nevertheless, although the global symmetry II is broken at the phase 
transition one can define a new operator n( s ) = exp (i-K^c + d^d]) which 
commutes with the superradiant Hamiltonian H s . 

The normal phase allows a ferromagnetic ordering for the pseudospins 
((J z ) — ► — j), whereas in the superradiant phase, (J z ) decreases continuously. 

1.3.4 Second-order quantum phase transition 

The excitations of the system are given by the energies E± , which describe 
collective modes or polaritons in solid-state physics. As the coupling A ap- 
proaches A c one can observe that = E^ = 0, signaling the occurrence of 
the quantum phase transition. In contrast, E + tends to the value ^Juj 2 + to 2 
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as A — ► A c from either direction. In fact, one may identify _E_ with the ex- 
citation energy of a photon branch and E + with the excitation energy of an 
atomic branch [51]. In addition, for A — > A c , from either direction one gets 
£7_(A — > A c ) cx |A C — A| z,y , where the dynamic critical exponent reads z = 2 
and v = 1/4 is the critical exponent describing the divergence of the charac- 

— 1/2 

tcristic length £ = E_ . The fact that vanishes at A c implies that the 
Dicke model exhibits a second-order quantum phase transition. The entan- 
glement between the atoms and field diverges with the same critical exponent 
as the characteristic length [53] . Emary and Brandes have also shown that at 
the quantum phase transition the system changes from being quasi-integrable 
to quantum chaotic [51]. The finite-size scaling exponents of the Dicke model 
have been discussed in Rcf . [54] . The quantum phase transition for the multi- 
mode (and continuum) case has been discussed in Ref. [55]. (Tolkunov and 
Solenov argue that adding the A 2 term in the Hamiltonian, one would observe 
the same transition with a corrected position of the critical point [55].) 

1.3.5 Realizations 

In familiar quantum-optical systems, the frequencies u and u s exceed the 
dipolc coupling strength by many orders of magnitude. Therefore, the quan- 
tum dissipation due to atomic spontaneous emission and cavity loss is usually 
unavoidable, and the quantum phase transition remains unobserved experi- 
mentally On the other hand, the Dicke model in the quantum phase transition 
regime {w s « lu ~ A) may be realized based on the collective interaction of an 
ensemble of atoms with laser fields and field modes of a high-finesse optical 
resonator. In particular, cavity quantum electrodynamics (QED) [56] might 
realize the Dicke model with parameters w A that are adjustable 

and can in principle exceed all dissipation rates [57]. More precisely, Dimcr 
et al. have proposed a well-defined scheme based on multilevel atoms and 
cavity-mediated Raman transitions to realize a Dicke model in an open sys- 
tem dynamics (with omission of the A 2 term) [57]. The ensemble of atoms 
is simultaneously coupled to the quantized field of the optical cavity mode 
and the classical field of a pair of lasers. In principle, optical light from the 
cavity carries signatures of the critical behavior. Another scheme based on a 
superconducting quantum interference device coupled to a high-quality cavity 
supporting a single-mode photon has also been proposed in Ref. [58]. 



1.4 Jaynes-Cummings lattice 



Finally, we consider a two-dimensional array of coupled optical cavities, which 
may be realized in circuit QED for example, each containing a single two-level 
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atom (spin) in the photon-blockade regime. The coupling between the atom 
and the photons leads to an effective photon-photon repulsion; this photonic 
repulsion (blockade) has been shown recently by Birnbaum et al. using a 
single trapped atom [59]. Below, we study the resulting Jaynes-Cummings 
[60] lattice-type model. Following a mean-field theory, first we show that the 
system at zero temperature can undergo a characteristic Mott insulator (exci- 
tations localized on each site) to superfluid (polaritons delocalized across the 
lattice) quantum phase transition [61]. Then, we make a rigorous comparison 
between the Jaynes-Cummings lattice and the Bose- Hubbard model [62]. 

1.4.1 Hamiltonian 

We consider a system composed of a regular array of identical cavities. For 
a sufficiently large quality factor, we may restrict our treatment to a single 
photon mode, similar to the Dicke model above. We assume that the finite 
quality factor only stems from photon leakage k among nearest-neighbor cav- 
ities. The Jaynes-Cummings lattice model then takes the following form [61]: 

H = E H i° - K E { a l a i + a ) a ) ~ ^E ( a l a * + a t^) > (1-36) 

i i 

where we have introduced the Jaynes-Cummings Hamiltonianll : 

H( c = Lua\ ai + w s o-fo-r + p(a\o-r + a + ai ). (1.37) 

The operator (aj) annihilates (creates) one photon in a given cavity at 
lattice site i. Similarly af denote the Pauli raising and lowering operators 
for each two-level atom. The bracket notation (i, j) denotes summation over 
nearest neighbor pairs. The phase boundary between a Mott insulator and a 
superfluid phase can be determined in a grand canonical approach. To this 
end, a chemical potential \x is introduced. The grand canonical approach con- 
siders a situation in which particle exchange with the surrounding is permitted 
and is used because of its convenience for determining the phase diagram. The 
boundary between Mott insulator and superfluid phases is determined by the 
value of [i for which adding or removing a particle does not require energy. In 
addition, here we work in the limit of the rotating-wave approximation. 

1.4.2 Mott insulator-superfluid transition 

In the limit n — > ("atomic" limit), the Hamiltonian decouples in the site 
index and reduces to a Jaynes-Cummings Hamiltonian H K=0 = H JC — /in 
where the total excitation number on each site now obeys n = a^a + a + o~~ . 
Since H JC and n commute we can write the eigen-energies of H K=0 as £\ n ±) = 



II The light-atom coupling f3 should not be confused with 1/ksT. 
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FIGURE 1.3 

(left) Ground state for k = and A = 0. (right) Mean-field phase diagram. 

E\n±) ~ A* n ) where the usual eigen-energies of the Jaynes-Cummings Hamil- 
tonian are given by: E\ n ±\ = nuj — A/2 ± y/ n(3 2 + A 2 /4 for n > 1 and 
i?| ) = for the ground state. Here, A = uj — iv s is the detuning between 
the two- level atom and the resonator frequency. The polariton states \n±) = 
(\n,i) ±|(n— 1), t))/\/2 are simultaneous eigenstates of the Jaynes-Cummings 
Hamiltonian and of the polariton number n = a* a + a + a~ . For a gien n, since 
the state \n+) is always higher in energy compared to \n—), we can ignore it 
completely when focussing on ground state properties. The transition from 
the boundary \n— ) to |(n+ 1)—) occurs when £\ n -) = £|(n+i)-) ( n — 0, 1, ••■)■ 
Therefore, this implies: (jji - w)//3 = y/ n + (A/2/3) 2 - v /n+T+ (A/2/3) 2 . 
This reduces to (fi — uj)/f3 — y/n — yn + l in the resonant case a; = w s . Those 
boundaries separate the different Mott phases for the polaritons (excitations) ; 
see Fig. 1.3 (left). These results for the atomic limit are consistent with the 
results of Ref. [61]. By increasing the hopping k, one expects a second-order 
quantum phase transition from a Mott-insulating (MI) phase to a superfluid 
phase (SF) of polaritons which can be first described using a mean-field theory. 

Introducing the order parameter ip = z c n{ai) where z c is the coordination 
number of the lattice, the phase boundary between the MI and SF phases can 
be determined in way similar to to the procedure for the Bose-Hubbard model 
[63]. In the critical region the ground state energy can be expanded as [62]: 

E G <S>) = Eg* + r#| 2 + iuH 4 + OM 6 ). (1.38) 
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This represents the standard situation of a quadratic plus quartic potential, 
ubiquitous in the study of mean-field phase transitions. At the mean-field 
level, the phase boundary is specified by the condition r = 0. In fact, the 
coefficient r can be expressed as r = R n + (z c k) _1 where R„ can be obtained 
from second-order perturbation theory in the photon hopping [61, 62] leading 
to Fig. 1.3 (right). Recently, a strong-coupling theory to the phase diagram 
has been developed in Ref. [64] which allows to include the leading correction 
due to quantum fluctuations for the phase boundary. The phase boundary is 
in agreement with Quantum Monte Carlo calculations in two dimensions [65] . 

1.4.3 Spin-1/2 mapping for the polaritons 

To make a more concrete analogy with the Bose-Hubbard model, one can build 
a spin-1/2 mapping for the polaritons between two adjacent Mott lobes. More 
precisely, the polariton operators are defined as [66] P} na — \na)i(0— |j, where 
a = ± and i denotes the site index. Note that these operators do not satisfy 
the canonical commutation relations of creation and annihilation operators. 
Now, close to a degeneracy point between two Mott-insulating lobes (polariton 
occupation numbers n — 1 and n) in the atomic limit k/ (3 <ti 1, one can build 
an effective model for low-energy states in this regime from the two relevant 
product states |(n— 1)— )® j and \n— )® J . Within this truncated Hilbert space, 
the Hamiltonian H^f? can be re-written as a spin-lattice XX model. For 
this purpose, one can introduce spin-1/2 operators cr 4 + = P} n _Pi ; („_i)_ [66]. 
Within this definition, one can check that erf and its Hermitian conjugate a~ 
obey standard commutation and anticommutation rules for the Pauli lowering 
and raising operators [62]. The XX model describing the physics close to the 
degeneracy points of the Jaynes-Cummings lattice model is given by [62, 66]: 

H* X = \ (£\ n -) -£|(»-i)->) 5>* (1-39) 

i 

The conversion amplitudes t n are given in Ref. [62] . This effective description 
is reminiscent of the Bose-Hubbard model between two Mott lobes. 

1.4.4 Field theory approach of the transition 

In fact, a field-theory formulation can also been built by analogy to the 
Bose-Hubbard model [62]. After performing a usual Hubbard-Stratonovich 
transformation to decouple the hopping term, then one obtains the following 
effective Lagrangian for the auxiliary fields ip*(x,r) and ip(x,T) (which are 
proportional to the order parameter (aj) and hence are small in the critical 
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region) : 



C eff =K + K,r^ + K 2 ^ +K 3 |V^| 2 + r\M 2 + Jvf- (1-40) 



dr 



dtp 



dr 



2 



As in the case of the Bosc-Hubbard model [67, 68] , the coefficients f and u can 
be related to the mean-field coefficients: f = w _1 r and u = v u where v is 
the volume per lattice site [62] . To obtain the coefficients K\ and K 2 one can 
use the fact that the theory must be invariant under a (time-dependent) U(l) 
symmetry. Using the gauge symmetry of the action S e ff = dr J d d xC e ff , 
then one obtains the precise (exact) equalities [62]: 

dr I <9 2 f 

K X = - and K 2 = --- -. (1.41) 

o(u) — fij 2 o(uj - fi) z 

Therefore, whenever the coefficient K\ vanishes, then the phase transition 
changes its universality class since the dynamical critical exponent z = 1. 
The physics of these multicritical points is in complete analogy to the cor- 
responding physics of the Bose-Hubbard model [63, 67, 68]. In contrast, for 
K\ 7^ 0, we have a rather different field theory and the dynamical critical 
exponent associated with the quantum phase transition is z = 2. 

On the other hand, the Jaynes-Cummings lattice has one additional pa- 
rameter, the energy scale for the atoms oj s . Thus, the phase boundary is a 
two-dimensional surface in the space spanned by the parameters u>, u> s , and 
k, and the condition K\ = defines curves on the phase boundary, whose 
position is completely determined by dR n j dui\ UJs = [62]. For each Mott 
lobe, there is one such multicritical curve, along which the universality class 
changes in a similar way to the Bose-Hubbard model. It is interesting to 
observe that the evidence for these multicritical curves is currently discussed 
controversially. Schmidt and Blatter [64] have presented evidence for the pres- 
ence of such multicritical points which is consistent with the results from field 
theory shown above [62]. However, Zhao et al. [65] argue for the absence of 
such multicritical points based on quantum Monte Carlo simulations. On the 
other hand, Zhao et al. have defined the dynamic critical exponent z from 
the superfluid stiffness of the photons only (and not of the polaritons); in 
particular, the photon number cannot be conserved (at the lobe-tips). 



1.4.5 Realizations 

Several candidates have been proposed for an actual realization of the Jaynes- 
Cummings lattice, ranging from arrays of photonic band-gap cavities to circuit 
QED systems [69, 70]. A natural candidate system is the microwave strip line 
resonator for circuit QED [71]. We consider that the realization in circuit 
QED is especially interesting [62]. Mostly, the basic building block is well es- 
tablished and show the required Jaynes-Cummings physics [72]. In particular, 
medium-size arrays (with a number of sites between 10 and 100) should not 
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pose fundamental difficulties and ideas to prepare the system have been envi- 
sioned [66]. First, the system would be prepared in the Mott insulating regime 
by using a global external microwave signal. Then, the system could cross the 
phase boundary by varying the detuning A, whereas the hopping strength k is 
fixed by the fabrication parameters. The realization of this Jaynes-Cummings 
lattice would share many of the fascinating aspects of the ultracold atoms. Fi- 
nally, a relevant question to address in the future is whether the effect of an 
additional external driving combined with the presence of dissipation may 
lead to a change of the universality class of the quantum phase transition. 



1.5 Conclusion 

In this chapter, we have provided a comprehensive and modern investiga- 
tion of quantum phase transitions emerging in spin-boson Hamiltonians. The 
thermodynamic limit here can be achieved via a large number of bosonic 
modes resulting in a dissipative environment or/and through a large ensem- 
ble of spins (spin array). In the case of a single two- level system coupled to a 
bosonic (dissipative) environment, we have summarized recent developments 
on the computation of observables and spin dynamics for the ohmic case and 
we have discussed the quantum phase transition in the sub-ohmic case by 
analogy to Ising spin chains with long-range interactions. Additionally, we 
have shown that a spin array coupled to the same dissipative bosonic bath via 
a longitudinal coupling gives rise to a dissipative quantum Ising model which 
may be engineered in cold atomic systems. On the other hand, the light- 
atom interaction also results in fascinating quantum phase transitions such as 
the zero-temperature superradiant transition taking place in the Dicke model 
when tuning the dipole coupling strength A. The Jaynes-Cummings lattice 
system comprising an array of optical cavities each containing a single atom 
also allows to realize a superfluid-Mott insulator transition for polaritons. The 
Dicke model and the Jaynes-Cummings lattice may be implemented in QED 
cavity systems. In this sense, setups of electromagnetic resonators could em- 
body a novel class of quantum simulators of condensed-matter Hamiltonians. 
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